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We present a method for detection of weak continuous signals from sources in binary systems via the in- 
coherent combination of many "short" coherently-analyzed segments. The main focus of the work is on the 
construction of a metric on the parameter space for such signals for use in matched-filter based searches. The 
metric is defined using a maximum likelihood detection statistic applied to a binary orbit phase model including 
eccentricity. We find that this metric can be accurately approximated by its diagonal form in the regime where 
the segment length is <K the orbital period. Hence correlations between parameters are effectively removed by 
the combination of many independent observation. We find that the ability to distinguish signal parameters 
is independent of the total semi-coherent observation span (for the semi-coherent span » the segment length) 
for all but the orbital angular frequency. Increased template density for this parameter scales linearly with the 
observation span. We also present two example search schemes. The first uses a re-parameterized phase model 
upon which we compute the metric on individual short coherently analyzed segments. The second assumes 
long » the orbital period segment lengths from which we again compute the coherent metric and find it to be 
approximately diagonal. In this latter case we also show that the semi-coherent metric is equal to the coherent 
metric. 



I. INTRODUCTION 



The search for continuously radiating sources in binary sys- 
tems has proven to be a consistently intensive endeavor in a 
number of branches of astrophysics. The prime focus being 
the search for pulsars in the fields of Radio and X-ray astron- 
omy, and more recently for non-axisymmetric rapidly rotat- 
ing neutron-stars in gravitational-wave (GW) astronomy. In 
the Radio and X-ray fields it has long been known how to 
perform searches for such systems over relatively short ob- 
servation times (<K orbital period), known as "acceleration" 
searches (T}(3). In this type of search the frequency of a 
Doppler modulated signal from a source in a binary system 
is, over the short observation, approximated by its Taylor ex- 
pansion and the search is performed over a parameter space 
defined by the frequency and its derivatives. More recent 
searches, sensitive in the complementary extreme to systems 
with orbital periods <K the observation time, have also been 
successful in the detection of pulsars in the Radio band 0- 
6|. These searches take advantage of the distinct frequency 
domain signature generated by a frequency modulated signal 
such as a Doppler modulated continuously emitting pulsar in 
a binary system. This type of search is also planned for appli- 
cation to the GW case Q. 

In this work we are motivated by the problems inherent 
to detection of GW radiation from known binary systems, 
of which primary examples are the low-mass X-ray binaries 
(LMXBs). In this case the signal is expected to be extremely 
weak and the parameter space known to be very large, the 
proverbial "needle in a haystack". This is also a long standing 
issue for X-ray astronomy for a subset of these objects [8] [9]. 
Here we consider known objects as systems where the sky- 



position is assumed to be known precisely and that the intrin- 
sic source frequency, as well as the orbital parameters of the 
system, are unknown. In the search for GWs from such ob- 
jects, to date, two strategies have been employed, one being 
a fully coherent analysis [ 1 1 spanning only a short observa- 
tion time, and the other employing a cross-correlation tech- 
nique [ 1 1 1 using data from multiple detectors. In the former 
a parameter space covering was used based on the coherent 
parameter space metric on a non-eccentric phase model with 
known orbital period l[T2l . The coherent metric, as we will 



show in the Section IC is simply a measure of "distance" de- 



'Electronic address: chris.messenger@astro.cf.ac.uk 



fined on the parameter space and informs us on how to place 
templates (or filters) optimally within the parameter space for 
a coherent search. Here we aim to expand on this approach 
and compute the semi-coherent metric, a similar measure of 
distance but defined on a semi-coherent detection statistic. 

The first work on the coherent parameter space metric for 
sources in binary systems [12] has recently been built upon 
by [ 13 1 where it has been shown for the GW case that searches 
for the LMXBs are computationally bound meaning that the 
number of templates one is required to process to optimally 
cover the parameter space is too large to computed on human 
time-scales. The apparently optimal fully-coherent approach 
is therefore unfeasible. This is a common theme within con- 
tinuous GW data analysis for large parameter space searches 
and has been addressed via various applications of what we 
will call "semi-coherent" searches. This is where an observa- 
tion is divided into a number of individual coherent observa- 
tions, or "segments", which are then incoherently combined 
to form a more powerful detection statistic. The term "semi- 
coherent" is in reference to the fact that the unknown phase of 
the signal has been either maximized or marginalized away 
within each segment of the analysis and hence there is no 
required phase coherence between the signal and our phase 
model for the duration of the entire observation span. Since 
the number of templates required for fully-coherent searches 
typically scale with the coherent observation time to a power 
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I > 1, by dividing up the observation into M segments the 
number of templates required to cover all segments (the num- 
ber for each segment multiplied by the number of segments) 
will be reduced by a factor of M ,_1 . The missing piece of the 
puzzle for sources in binary systems is the method by which 
to combine the detection statistics computed on the templates 
within each coherently analyzed segment. This is also true in 
the case of the search for X-ray pulsations from the LMXBs 
where, in addition to being computationally bound just as in 
the GW case, single observations of individual sources from 
X-ray timing satellites such as RXTE (Rossi X-Ray Timing 
Explorer) are typically restricted to time-scales of order ~ 1 
hour. This limitation, due to scheduling constraints and obscu- 
ration of the source as the detector orbits the earth, represents 
an automatic and unavoidable division of the complete dataset 
into segments. 

The work presented here is based on that of lfT4l4T6l who 
were the first to describe the concept of the metric, based on a 
distance measure equivalent to the expected loss in the detec- 
tion statistic between signal and template (the mismatch) in a 
matched-filter based search. From such a metric one would 
then proceed to place templates within the parameter space 
based on the criteria that the mismatch at any given point was 
always less than a desired threshold. In ifTTl is was shown 
that the metric on the general semi-coherent detection statistic 
was simply the average of the coherent metrics from the con- 
stituent segments. Recent work by |[T8l[T9ll has described the 
construction of such a metric on the parameter space defined 
by the frequency and it's derivatives and the sky coordinates. 
This result is specific to all-sky, wide band, searches for GWs 
from unknown isolated sources. In this work we compute the 
semi-coherent metric for the complementary class of sources, 
the known sources with known sky position in binary systems 
with unknown frequency and orbital parameters. 

In the remainder of this section we review the binary sys- 
tem signal model and the power-like detection statistic as well 
as the coherent and semi-coherent metric definitions. In Sec- 
tion [n] we describe the calculation of the coherent and semi- 
coherent metric for the specific case of sources in binary sys- 
tems. In particular we analyze two complementary scenarios, 



in Section II A we deal with cases where the segment length is 



<K the orbital period of the source and in Section II C the case 



where the segment length » the orbital period. In each case 
we outline a basic method for the practical implementation of 



the intrinsic source frequency is constant (i.e. there is no in- 
trinsic frequency-evolution of the source) and that any detec- 
tor motion relative to the chosen inertial reference frame has 
been accounted for. The latter assumption in practice requires 
that the sky-position of the source is known to high accuracy 
and that the time-series has been barycentered (usually to the 
solar-system barycenter). We also ignore the slowly varying 
amplitude response in the case of a GW signal. The remain- 
ing phase contribution is assumed to be entirely due to the 
constant intrinsic frequency and the Roemer delay across the 
source orbit. 

This time-dependent phase for a source in a bound eccentric 
(non-relativistic) orbit is |20| 

O(f) = 27rv{f-f ref (2) 
+ a [sin w(cos E(t) - e) + cos a> sin E(t) Vl - e 2 J 

where v is the intrinsic signal frequency 1 , f re f is some refer- 
ence time (at which the signal phase is equal to <E>o in the 
source frame), a is the orbital semi-major axis projected along 
the line of sight and divided by the speed of light, u> is the ar- 
gument of periapse, e is the orbital eccentricity and E is the 
eccentric anomaly. The eccentric anomaly is defined by the 
transcendental relation 



In 



— (t - f p ) = E - e sin E 



(3) 



where P is the orbital period (equal to 2n/Q. where £2 is the 
orbital angular frequency) and f p is the time of periapsis (the 
point of closest approach to the prime focus of the orbit). 

With regards to the search for continuous GW emission 
from sources in binary systems, the primary targets (as men- 
tioned in Section [I) are the LMXBs. These sources are ob- 
served to have highly circularized orbits (see iTPSIl for de- 
tailed descriptions of current orbital parameter estimates) and 
as such we have chosen to limit our investigations to low- 
eccentricity orbits. In the low-eccentricity limit (e <K 1) we 
can adopt the approximation BT1I221 



<D(f) * 2nv\t-t mi 



(4) 



a search, and in Section III we summarize our findings. 



+ a 
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sin ifr(t) + - sin 2ifr(i) - + cos 2tff(t) - — 



A. The signal model 

We model a continuously emitting source of radiation (ei- 
ther GW or electromagnetic (EM)) located within a binary 
system with a non-emitting companion. In this sense we de- 
fine the noise-free continuous signal received at an inertial ref- 
erence frame to be 



s(f) = A sin [<D(r) + O ] 



(1) 



where A as the constant signal amplitude and Oo is a constant 
phase offset. Note that in the following we are assuming that 



where we have Taylor-expanded the orbital contribution to the 
phase in powers of the eccentricity e up to leading order. In 
addition we have adopted the following parameters, 



k — e cos a>, 
j] = esinw, 



(5a) 
(5b) 



1 We have assumed that the velocity of the binary barycenter is constant 
relative to the inertial reference frame and hence absorbed any Doppler 
shifts into v. 
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as replacements for the more physical e and a> parameters. 
This change is motivated by the fact that for low-eccentric 
systems strong degeneracies between the argument of peri- 
apse co and the time of periapse passage f p would complicate 
our analysis. In addition we also define 



(6) 



as the time dependent orbital phase as measured relative to 
the time of passage through the ascending node of the orbit 
f asc = t p - to/Q.. 

The circular orbit case is a specific instance of the more 
general elliptic orbit (Eq.|2]). For e — we can write the phase 
model in this case as 



<b{t) = 2nv [t - f ref - a sin (//(?)] 



(7) 



where the extraneous parameter a> has been set to to — n such 
that, just as for the Taylor-expanded low eccentricity case, the 
orbital reference time is the time of passage through the as- 
cending node of the orbit rather than the, now meaningless, 
time of periapsis f p . 

For the majority of our primary sources, the LMXBs, the 
expected eccentricities in these systems are relatively low 
(typically > 10~ 3 ). However, despite the apparent current re- 
quirement for circular orbits only, we continue the analysis 
using the phase defined for an low-eccentricity orbit (Eq. |4]i 
since the results that follow from this choice can easily be in- 
terpreted for either the eccentric or circular orbit cases. 



B. The detection statistic 

In searching for deterministic signals in noisy data one 
wishes, in general, to distinguish between the noise hypothe- 
sis and the signal-plus-noise hypothesis in order to determine 
the presence of a signal. In this case the likelihood-ratio test is 
optimal in the Neyman-Pearson sense 11231 . If we assume that 
our measured time-series data-set consists of our signal and 
additive Gaussian noise sampled uniformly at discrete times 
tj such that 



x(tj) = s(tj, 0') + n(tj), 



(8) 



where 0' is a vector containing the signal parameters. We can 
write the likelihood function as 



1 



N N 



l(0') <* ex P \ - - 2 2 ( Xj - sjmy-ji ( Xk - Sk m) 



(9) 



where C is the noise covariance matrix, N is the number of 
samples in our time-series and for simplicity of notation we 
have replaced x(tj) with xj and s(tj) with sj. 

One can show that twice 2 the log-likelihood ratio (the log 
of the ratio between the likelihood function defined above and 



the likelihood assuming no signal), when analytically maxi- 
mized over the "nuisance" parameters A and <t>o, becomes 



A(0) = 



S n AT 



N 



(10) 



in the case of "white" noise where the covariance matrix is 
diagonal 3 . Note that the un-primed represents the reduced 
parameter set after having removed the dependence upon the 
nuisance parameters and that we have used At as our sampling 
time, AT = NAt as our time span and S n as the single-sided 
noise spectral density (assumed constant over the frequency 
band of interest). This detection statistic has the form of a 
Fourier power with the exception that the complex phase, with 
which we multiply each datum, contains an orbital phase com- 
ponent in addition to the intrinsic frequency component. 

The statistical behavior of the log-likelihood ratio detection 
statistic in additive Gaussian noise follows that of a x 1 dis- 
tribution. In the simplified constant amplitude case described 
here the random variable A(0) is described by a non-central 
x\ distribution (the subscript indicates the number of degrees 
of freedom). Assuming a set of signal parameters the ex- 
pectation value of A(0) evaluated at an offset parameter space 
location + AG is then given by 



E[A(0,A0)] = 2+p 2 (0,O) 



Z 



(ID 



where p 2 (8, 0) is the optimal signal-to-noise ratio (SNR) given 
by 



P 2 (0,0) 



Sn 



|>4(0). 



(12) 



7=1 



Note that we have used A<£ ; (0, AO) = <J> ; -(0+A0)-O/(0) to rep- 
resent the phase offset caused by the offset in parameter space 
location. Also note that the second term in Eq. [TT] is equal 
to the non-centrality parameter governing the j^distribu- 
tion which becomes the optimal signal-to-noise-ratio (SNR) 
squared when the template and signal parameters are exactly 
matched. In our problem we are faced with a continuum of 
signal hypotheses defined by the parameter space spanned 
by 0, which we must somehow sample from. The common 
frequentist strategy in parameter space searches is to maxi- 
mize the detection statistic over all unknown parameter val- 
ues. For two of these parameters, A and Oq, we have per- 
formed this maximization analytically. In the following sec- 
tion we describe the standard strategy for performing the re- 
maining maximizations 4 . 



2 We multiply the log-likelihood ratio by 2 so that it is exactly consistent 
with ax 2 distribution with 2 degrees of freedom. 



3 This is a valid assumption for practical purposes since analyses are typi- 
cally divided into narrow frequency regions within which the noise spectral 
density can be assumed "white". 

4 We note that strictly speaking, the optimality of the standard likelihood- 
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For a deterministically amplitude modulated signal, such as 
in the GW case, the standard technique is to analytically max- 
imize not only over the phase and amplitude of the signal but 
also over the inclination angle of the source and the polariza- 
tion angle of the GW wave (25 1. In such a case the detection 
statistic becomes a^ 2 statistic with 4 degrees of freedom. 



C. The coherent metric 

The parameter space in our case is defined by the ranges in 
uncertainty on the frequency and orbital parameters defining 
our signal. When faced with the prospect of searching this 
space for the true signal parameters we rely on the concept of 
the parameter space metric lfT4l - [T6ll which allows us to set a 
measure of distance by which we can determine how to sam- 
ple within the space. Using this geometrical approach we are 
able to satisfy the constraints that our templates will not be 
placed too coarsely such we will not "miss" a signal and also 
that they will not be placed too finely such that we will be 
wasting computational effort. The standard choice is to de- 
fine the distance measure for a coherent analysis as the ratio 
of the expectation value of the loss in SNR local to a signal's 
true parameters and the expectation value of SNR at the true 
signal parameters. This measure, or mismatch, is then 



(i(0, A0) = 



p 2 (e,O)- P 2 (o,A0) 

P 2 (0,0) 



(13) 



where we use By Taylor-expanding the mismatch around the 
true signal location we obtain 



H{6, AO) = g^QSf) Aff"A8 v + o(A6> 3 ) 
where the coherent metric g^ifl) is defined as 

fm _ l m«) wo) \ l \ l wfl \ 

WBJ "\ 80" d& I \ dp l\ 36* I 



(15) 



with (...) representing the average over the coherent observa- 
tion time and with A0^ representing the deviation in the fi'th 
parameter from its true value. 



D. The semi-coherent metric 

The first stage of a semi-coherent analysis is the process of 
performing the multiple constituent coherent analyses on the 
M independent data segments into which the full observation 
has been divided 5 . Such separate analyses will result in the 



ratio applies only for the point-hypothesis case. It has been shown 1241 
that the optimal statistic in the more general case (including the point hy- 
pothesis case) where one is faced with a continuum of signal hypotheses 
(defined one some parameter space) is the Bayesian statistic, the "Bayes 
Factor". 

5 Whilst in this work we concentrate on a division of the complete dataset in 
the time domain, semi-coherent searches in general may also suit division 
of the dataset in the frequency domain, e.g. 1 7 j . 



generation of discretely sampled values of the detection statis- 
tic A,„ on the parameter space spanned by 0, where m indexes 
the segment number. Note that the metric in each segment and 
therefore also the mismatch will vary between segments since 
the metric can be a function of the segment epoch (as is the 
case in our binary system). 

Since A,„ are maximized log-likelihood ratios it follows 
that a sensible choice is to define the semi -coherent detection 
statistic as the sum of these values as a function of 0, 



M 

A(0) = ^A m (0). 

m=l 



(16) 



Therefore, such a statistic would itself represent a log- 
likelihood on the space spanned by 0, with the understanding 
that there has been an implicit maximization over M distinct 
amplitudes A m and d>o,»i- 

In practice, the reason behind adopting an semi-coherent 
strategy is likely that we lack the computational resources 
required to coherently track the phase of a signal over long 
observation times. In which case the introduction and max- 
imization of M distinct initial phases is an approximation to 
the true phase model. We note that the similar introduction 
of the M distinct signal amplitudes also does not appear to be 
consistent with our original signal model. The semi-coherent 
detection statistic is now also sensitive to signals for which 
the amplitude varies with a time-scale » the coherent segment 
length AT. This subtle feature is a potential improvement with 
regards to searches in the X-ray spectrum for signals from 
LMXBs where the amplitude is not necessarily expected to 
be constant [26 1 . 

Based on the fact that the individual segment A m values are 



(14) x\ distributed it follows from Eq. 



16 



that A is x\ M distributed 
with a non-centrality parameter equal to the sum of the in- 
dividual segment non-centrality parameters. The expectation 
value of A(0) is therefore given by 



M 

E[K(0,A0)] =2M + 2p 2 (0,A0). 



(17) 



The semi-coherent mismatch (defined as the loss in semi- 
coherently summed SNR) is then 



(40, AO) 



1 M 

-T 



m=l 
M 



p 2 m (0,0)-pl(0,A0) 

rin(0,o) 



(18a) 



(18b) 



m=l 



where we have assumed that each segment has identical du- 
ration and noise level. In practice, whilst it may be simple to 
maintain constant segment lengths, detector noise may vary 
between segments. In addition, for GW detectors, amplitude 
modulation of the signal due to the changing response of the 
detector as the earth rotates will have an effect ll27l . 

By substituting Eq.[l4]into Eq. |18b| the semi-coherent met- 
ric Gfj V (0) can then be expressed as 



Guy(0) 



1 M 



(19) 
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This is the standard result, as shown in [17|, that the semi- 
coherent metric is simply the element-by-element average of 
the constituent coherent metrics. 



E. Template placement 

In both the coherent and semi-coherent stages of the anal- 
ysis the intention is to use the information contained in the 
metric description of the parameter space to place templates. 
The task is then to cover such a space efficiently whilst adher- 
ing to the constraint that no template should have mismatch 
with any potential signal, greater than a given threshold, fi*. 
In this work we simply refer the reader to recent efforts made 
in the field of GW data analysis with regards to metric based 
template placement 6 . 

In the case of a constant metric, one in which either care- 
ful parameterization or luck has left a metric of which the el- 
ements are constant over the parameter space, the problem 
is directly equivalent to the "sphere covering" problem (28). 
This was realized and investigated, in the context of GW data 
analysis, in (39). In general, in the constant metric case the 
optimal solution (in the sense of guaranteed parameter space 
coverage) is to use an n-dimensional lattice of templates. The 
most basic (but inefficient) of these lattices being the hyper- 
cubic Z„ lattice and the most efficient being the class known 
as the A* lattice. 

For the more general case of a metric whose elements are 
functions of the location in the space itself, template place- 
ment is more difficult. The coordinate volume and orienta- 
tion of the area local to a template will change as one moves 
through the parameter space. Hence the density and relative 
separations between templates in each of the dimensions of 
the space changes depending on where we are placing tem- 
plates. An optimal solution to this particular problem has yet 
to be found although recent work on so-called "random" and 
"stochastic" template banks O0H331 shows impressive cov- 
ering performance especially for higher dimensional spaces. 
These methods use only information regarding the required 
local template density, a quantity proportional to the square 
root of the determinant of the metric. 



far more finely than the coherent template banks, and b) poten- 
tially defined on a different parameter space coordinate system 
(this is the case for the scenario described in Section |lTA"] >. For 
each semi-coherent template we need to perform the follow- 
ing. 

Firstly, for each coherent segment we transform the param- 
eters of the semi-coherent template into the parameters used 
in the segment. For each set of transformed coordinates we 
interpolate the value of the detection statistic on the coherent 
template bank at the desired parameter space location. The 
interpolated detection statistics are then summed according to 



Eq. 16 This is repeated for each template in the semi-coherent 
template bank. 

A particularly fast and simple interpolation method is to 
use the nearest neighbor approach, the speed and simplicity 
of which benefits from the use of a hyper-cubic lattice within 
each segment. Higher precision is achieved by using smaller 
mismatches in the coherent stage. Using both the inefficient 
hyper-cubic lattice and decreasing the mismatch will of course 
increase the computational cost of computing the coherent 
stage of the analysis. 

The assumption that the coherent analysis had been com- 
puted prior to this procedure need not be the case in prac- 
tice. One could perform the interpolation and summing as 
each segment is being analyzed as long as the total time span 
of the complete data-set was known a priori, since as we will 
see, the semi-coherent metric is dependent upon the total time 
span. 



H. RESULTS 

In this section we report on the results obtained through cal- 

i using 



culation of the semi-coherent metric elements via Eq. 19 



the binary system phase model given by Eq. [4] In parallel we 
provide basic schemes for the practical application of a semi- 
coherent search in two observational regimes determined by 
the length of the coherent observation AT. The first having 
a "short" coherent stage with AT <K P where the segment 
length is far shorter than the orbital period, and the second 
having a "long" coherent stage AT » P where the segment 
length exceeds the orbital period. 



E Combining results from different segments 

Let us imagine that the coherent detection statistic has been 
computed for all templates in all segments. We now have 
the information required to compute an ensemble of semi- 
coherent detection statistics and the general form of the semi- 



coherent metric (Eq. 19 1 tells us how to construct a semi- 
coherent template bank on the same parameter space. We note 
that the semi-coherent template bank is a) potentially sampled 



' Note that there are many other parameter space exploration strategies that 
do not use the metric e.g. Markov-Chain-Monte-Carlo methods. 



A. The coherent metric for "short" (AT <k P) coherent 
segments 

In practice one need not compute the coherent detection 
statistic for all segments using templates given by the semi- 
coherent metric. The semi-coherent metric tells us how finely 
we must sample the parameter space for the complete semi- 
coherent analysis. It is the coherent metric that tells us how 
finely we must sample each segment. 

In the "short" segment regime the coherent metric, calcu- 
lated in our physical frequency and orbital parameter space, 
develops strong parameter degeneracies for short observation 
times, making template placement difficult. In this situation 
the reduction in detection statistic due to offsets between the 
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template and the signal in one parameter can be effectively 
counteracted by offsets in another. Additionally, the coherent 
metric computed in any chosen coordinates, would in gen- 
eral result in different physical parameter template locations 
in each segment. Therefore, once the coherent segments have 
been processed, for a given semi-coherent template there will 
not be a collection of corresponding coherent segment detec- 
tion statistics to sum together all computed at the exact same 
parameter space location. 

To address the first issue we propose the adoption of a re- 
parameterization 1 1 8 1 that greatly simplifies the phase model 
but that is only valid in the limit OAr <K 1, i.e. when the 
length of coherent segment is only a fraction of the orbital pe- 
riod. In this case, by simple Taylor-expansion of the phase 
model, given by Eq. |4j in terms of the time t about the mid- 
point of each coherent observation we obtain an approximate 
phase model given by 



<&Jtj,u) 



k=i 



(m) 
±_ 

k\ 



(0 



. f (m)V 
mid/ 



(20) 



where t ™} A is the midpoint of the m'th segment, <t>l m) is the 



'mid ^ """F"" 11 >^ 6 "— •*"() 

phase at this midpoint and the new coordinates uT 1 map to the 
physical parameters via the equations given in Appendix [A] 
Note that the u coordinates themselves are exactly the instan- 
taneous phase derivatives at the midpoint of each segment in- 
dexed by k which runs from 1 to n. Also note that by this 
re-parameterization the boundaries of the parameter space in 
the new coordinates will not have the same shape as the phys- 
ical parameter boundaries. 
By direct application of Eq. 



15 to our re-parameterized 



phase model we find that the metric in the u coordinates is 



7T 2 pvAT^ +V 



g M v(u) = \ l^Hp + l)!(v + l)!(yU + V + 1)' 
0, 



p. + v - even 



p + v — odd. 
(21) 

Here we see that the elements themselves are independent of 
parameter space location making this a constant metric and 
therefore simple for template placement. The fact that the 
gfj V ( u ) elements vanish for p + v being odd indicates a lack of 
correlation between adjacent phase expansion terms, a feature 
that is exemplified in Fig.[T] In this figure we show the results 
of a simulation in which the coherent mismatch p(u, Ah) has 
been computed in a region in the u parameter space surround- 
ing a simulated noise-free signal. The degree of correlation 
between the parameters for the cases where p + v is even is 
clearly indicated by the diagonal orientation of the mismatch 
contours. Also evident from the figure is that the predicted 
10% mismatch contour, computed using the metric approxi- 



mation given in Eq. 21 is in good agreement with the simu- 
lation results. The metric itself is a quadratic approximation 
around the peak of the log-likelihood ratio and in general is 
only expected to be accurate in the regime p <K 1. As we 
move away from the true signal location the simulated con- 
tours begin to deviate from their elliptical shape as higher or- 
der contributions to the mismatch become important. Fortu- 
nately, since by design, the mismatch is the loss in expected 



SNR one would typically only ever place templates according 
to mismatches (9(10%) so as to recover a large enough frac- 
tion of the SNR to avoid missing any signals. We note that in 
this region the predicted mismatch using the metric approxi- 
mation would have discrepancies relative to the true mismatch 
at a level <0.1%. 

The re-parameterization of the phase into the u coordinates 
is in fact entirely equivalent to the approximations made in so- 
called acceleration searches OKU- Here, in just the same way 
as in an acceleration search, the length of coherent observa- 
tion coupled with the orbital period determines how many or- 
ders of expansion are required to accurately model the phase. 
Conversely, given a computationally limited number of expan- 
sion terms, there is a corresponding limit to the shortest orbital 
period that can be searched for a given coherent observation 
time. Based on worst case values of f asc in relation to we 
find that 



(OAF) < 



A<D(« + 1)! 
2nav 



(22) 



represents a limit on the length of the coherent observations 
for an n'th order expansion in the u coordinates where AO 
is the allowed error in phase between signal and model. For 
example, in order to model a v = 100 Hz signal with a pro- 
jected semi-major axis a — 1 second with error in signal phase 
A<X> = 7r/2 using only an n — 2 expansion we would be limited 
to coherent observations spanning up to *l/25 of an orbit. 

The number of templates required to cover a parameter 
space is proportional to its proper volume, the integral of 
the square-root of the metric determinant over the parameter 
space, and is given by 



N = t(n,p) 



(23) 



where we have used S to represent the parameter space and 
the template density %(n,p) is a function of the dimensionality 
of the space n, the desired mismatch p and the choice of lattice 
covering. This gives us 



«(«+l)/2 



d n u 



(24) 



as the number of templates required for a single segment. 
Note that we do not give the explicit scaling associated with 
the parameter space boundaries (as done below in Eqs. |28|29| . 
The integral over the parameter space volume is a complicated 
function of these boundaries and will behave differently de- 
pending upon the relative sizes of the uncertainties in the or- 
bital parameters e.g. well known orbital period with poorly 
known time of ascension compared to poorly known period 
and well known time of ascension. One of the simpler cases 
one can consider is one in which the time of ascension is com- 
pletely unknown. In this case the parameter space boundaries 
in the u space form a hypercube with limits v m i n < u\ < v max 

and -flmaxVmax^iax < u k < flmaxVmax^max for k > 1 and the 

number of templates is given by 



N™Kt{n,p)AT 



n(«+l)/2 



(Vmax - Vmin) (« 



max y max 



yi-1 «(«+l)/2-l 
) ^ ^max 

(25) 
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FIG. 1 : Shown here are the results of a simulation in which the mis- 
match fi(u) has been computed over a region in the 3-dimensional 
parameter around the location of a simulated signal. The simulation 
models a single segment of 200 seconds length with signal parame- 
ters v = 200 Hz, a = 1 second, fl = 10~ 4 rads s , e = 0.02 and 
a> = 1 rad (r asc chosen randomly from the range [-P/2, P/2]) Each 
panel shows the mismatch as a function of each pair of parameters 
where the mismatch has been maximized over all other parameters. 
The solid black and grey contours indicate the measured mismatch 
ranging from /x = 0.1 (black) in steps of 0.1. The dashed red el- 
lipses are those calculated using the metric given in Eq. |2T] on the 
approximate phase given by Eq.|20| 



Here we would like to highlight a general feature of the re- 
parameterization process. In simplifying the template place- 
ment procedure, in the sense of using a constant metric, we 
have introduced parameter space boundaries in the u space 
that both vary between segments and are in general no longer 
simple hyper-cubic spaces. 



B. The semi-coherent metric for short (AT <<c P) coherent 
segments 

In order to compute the semi-coherent metric in this case 
we make the following assumptions and choices. We assume 
that we have many individual "short" coherent observations 
of duration AT randomly and uniformly distributed over the 
total observation span t, and that t » P such that the total 
span of all observations contains many source orbits. We also 
assume that the time of passage through the ascending node 
of the orbit f asc has been chosen such that it is close to t — 
(defined as the midpoint of the span of all observations). This 
is always possible since f asc can be redefined for any orbital 
cycle by adding or subtracting integer multiples of the orbital 
period 1 . 

The resulting semi-coherent metric elements are given in 
Appendix [C] for the case of e «: 1, r » AT and t » P. We 
also give the approximate Taylor-expanded versions of the el- 
ements applicable in the limit that Q.AT <K 1 where the seg- 
ment lengths are far smaller than the orbital period. It is on 
these Taylor-expanded results that we now focus on. Note 
that throughout, Taylor expanded metric elements are labeled 
with the superscript (T). 

The specific parameter vector for which our metric has been 
calculated is 6 = {v, a, f asc , £2, k, if] and in order to compute the 
semi-coherent metric itself we use the following approxima- 
tion 

1 M 

k=l 

l r T/2 

* lim - £ v \e,t mid )dt mid (26) 

"™ T J-t/2 

where we have replaced the sum over segments with an in- 
tegral over the midpoint of each segment (approximated as a 
continuous variable in the regime where M » 1). Under this 
approximation the semi-coherent metric can be accurately ap- 
proximated as 



-(T) 
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(27) 



Note that we have made the approximation that all off- 
diagonal elements can be set to zero, indicating that there are 
no parameter correlations. This approximation can be vali- 
dated by inspecting each element pair in turn and noting that 



the off-diagonal correlation terms are, in all cases, negligi- 
ble. A quantitative validation of this statement can be per- 
formed by computing the volume element (the square-root of 
the metric) on each 2-dimensional subspace and noting that 
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the inclusion of the off-diagonal terms has insignificant effect. 
Significant correlations would reduce the volume element and 
ultimately lead to a reduction in the number of templates be- 
cause individual templates would occupy more space. If we 
consider, for example, the parameters v and a we find that 
by including the off-diagonal term in the calculation of the 
sub-volume element the fractional reduction in the volume is 
xa 2 Q. 2 /4 which is negligible for all practical orbital parameter 
values. 

The reason that the metric becomes approximately diago- 
nal and that therefore parameter correlations are removed can 
be explained as follows. The analysis of each segment can be 
viewed as an independent measurement of the signal and its 
parameters. Recall that we use templates placed in the u space 
rather than the physical parameter space 9 because of degen- 
eracies in the physical parameters for short observation times. 
This simply means that our ability to determine certain phys- 
ical parameters and certain physical parameter combinations 
is poor in this choice of coordinates. 

For example, imagine a segment coinciding with a small 
fraction of the orbit that places the source at the ascending 
node of the orbit. In this case the source is moving towards us 
at its greatest speed. We will also focus on a single parameter 
pair, the intrinsic frequency v and the projected orbital semi- 
major axis a but note that the proceeding arguments apply 
equally well to all other parameters. At this segment epoch 
we are essentially measuring the maximally Doppler modu- 
lated signal frequency a;v(l + aQ) Hz and we therefore find 
that there exists a degeneracy making us unable to disentan- 
gle the true value of the frequency from that of the projected 
semi-major axis. The signal can equally well be estimated 
within this segment as having a higher frequency and smaller 
semi-major axis or vice versa. In this sense we would say that 
there exists a negative correlation between these parameters. 
A negative (or positive) correlation implies that an increase in 
mismatch from a positive offset in one of the parameters can 
be counteracted by a negative (positive) offset in the other. 

If we now examine the behavior of the v, a correlations at 
an orbital epoch advanced in time by P/4 seconds we find that 
the source is now moving perpendicular to the observer's line 
of sight. In this case the frequency is approximately constant 
for the duration of the segment and there would be little infor- 
mation in the signal regarding it's orbital acceleration. There 
now exists no correlation between v and a and one is able to 
determine the frequency quite accurately at the expense of rel- 
atively poor determination of the projected semi-major axis. 
Advance again by P/4 and the source is now receding from 
the observer at its maximum speed. We now see exactly op- 
posite behavior to that seen for the initial epoch in that there 
now exists an equal in magnitude positive correlation between 
v and a. Finally, advancing again by P/4 we find that the cor- 
relations have vanished just as for the epoch one half orbit 
earlier. 

Under our assumptions that the total observation spans 
many orbits and that the short observations themselves are 
uniformly distributed over this span we obtain measurements 
from all orbital phases. The combination of a pair of measure- 
ments that have approximately equal and opposite correlations 



will therefore remove those correlations. With many segments 
randomly positioned in orbital phase each one can be approxi- 
mated as having a corresponding segment on the opposite side 
of the orbit. As such, each opposing pair of segments removes 
the correlations in the final semi-coherent measurement. 

Shown in Fig. [2] are the results of numerical simulations 
whereby the mismatch fi(0, AO) has been computed over the 6- 
dimensional parameter space surrounding a simulated signal. 
Overlayed on the plots within this figure are the analytically 
derived mismatch ellipses corresponding to solving Eq.[T9]for 
fi — 0. 1 using Eq. 27 as the semi-coherent metric. The simu- 
lation comprised of M — 500 coherently analyses AT = 200 
second long segments randomly distributed across an full ob- 
servation span t = 10 8 seconds. The simulated signal has an 
orbital angular frequency Q = 10~ 4 rads s _1 (orbital period of 
17.45 hours), a frequency v = 200 Hz, an orbital semi-major 
axis a = 1 second and an eccentricity e = 0.02 and argument 
of periapse a> - 1. It is clear that in this case the approxi- 
mations made in the Taylor expansion and discarding of the 
off-diagonal metric elements were valid since there is no indi- 
cation of correlation between parameters. Correlations would 
appear as ellipsoidal mismatch contours with semi-minor and 
semi-major axes inclined with respect to the parameter axes, 
i.e. tilted ellipses as seen in some panels of Fig.[TJ We also see 
again at large fi > 0.1 mismatches, effects due to the break- 
down of the quadratic approximation used in the definition of 
the metric. These deviations from ellipticity, seen the simu- 
lated mismatch contours, can also be attributed to the finite 
number (M = 500) of segments randomly distributed over the 
observation span. This is due to the fact that the mechanism 
through which one can assume the disappearance of the off- 
diagonal metric elements is applicable only in the limit where 
the number of segments M » 1. We note however, that in 
the region of interest around fi ~ 10% the discrepancy be- 
tween the metric approximation and our simulation is <1%. 
To place this in a practical context, using this metric to place 
templates with an expected maximum mismatch between sig- 
nal and template of 10%, in a worst case scenario could result 
in a true mismatch of *1 1%. 

Using the metric result we can compute the number of 
semi-coherent templates required for a search and the scalings 
associated with the parameter space dimensions. Considering 
the full n — 6 low-eccentricity orbit system we obtain 



7Vsemi = i6 Voy at6t ^ ~ ^ ( a ™ x ~ """" 



) 



X (^ max ^min) (^asc.max ^asc.min) ^maxi (28) 

where we have defined the parameter space S as a 6- 
dimensional cube bounded by the minimum and maximum 
ranges of the parameter vector 0. We have used the relation 
dKdrj - e\ cos(2aj)\d edai in computing the integral over the ec- 
centricity parameters. We have used lower and upper bounds 
of zero and e max respectively for eccentricity and assumed 
that the argument of periapse would always be completely un- 
known and hence would have a range of zero to 2n. Note that 
for the template density we have assumed sub-optimal tem- 
plate placement based on a hyper-cubic lattice. A reduction 
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of a factor of w 6.8 can be achieved if an A* lattice is used in 
this n — 6 case. 

We show in Fig.|3]the regions of parameter space for which 
a semi-coherent search would require coverage in the eccen- 
tricity parameters k and tj. This requirement is strongly depen- 
dent on the length of the constituent coherent observations, 
the product of the frequency and semi-major axis, and the ec- 
centricity itself. As stated earlier, for the majority of our pri- 
mary sources, the LMXBs, the expected eccentricities in these 
systems are relatively low (typically <10~ 3 ). In addition, the 
computational limitations imposed by the size of the remain- 
ing parameter space forces the coherent observation length to 
shorter times. This leads us to the conclusion that in general 
a circular orbit phase model will be applicable to the LMXB 
sources. For a circular orbit system, ignoring the k and rj di- 
mensions, we have the following template scalings 

NZu = ^^ Ar4r ( v -- v ™)K-- fl -) 

* (^max — ^min) (/asc.max — fasc.min) • (29) 

Note that again we have assumed a sub-optimal hyper-cubic 
template lattice and in this n — 4 case a factor of m 2.8 re- 
duction in templates can be achieved using an A* lattice. We 
wish to stress that common mistakes made in estimating the 
required number of templates stem from the inclusion of di- 
mensions that are geometrically "thinner" than the mismatch 
coverage of a single template. Such mistakes can lead to se- 
vere underestimation of the number of required templates. For 
example, take the parameters of the simulations for which re- 
sults are shown in Fig.[2]and apply them to Fig.[5]which shows 
the requirements of a semi-coherent search with regards to 
placing templates in the eccentricity parameters . We see that 
for Q.AT = 0.02 and va = 200 the corresponding eccentricity 
threshold is » 0.1. We would therefore expect a semi-coherent 
search t be unable to distinguish eccentric orbits from circu- 
lar orbits for eccentricities below this value. Looking at the 
mismatch ellipses shown in Fig. [2] we see that this is indeed 
the case, the scale of the ellipses in the k and 77 dimensions 



(k 2 + rj 1 - e 2 ) is ~ 0.1 whereas the signal has eccentricity 
e = 0.02. This is a clear example of when one should as- 
sume a circular orbit model and remove the eccentricity and 
argument of periapse from the analysis. 

The linear dependence of the number of templates on the 
total observation span t in both the eccentric and circular obit 
cases indicates that an increase in observation span results in 
an increase in the number of templates. One can view this 
as a "refinement" of the semi-coherent template bank with in- 
creasing observation span. For our case this refinement comes 
from a single metric element, Gnn, and governs the template 
spacing in orbital angular frequency alone. The reason that 
we see this behavior is that the way in which we combine seg- 
ments in this regime is a function of the orbital angular fre- 
quency (or equivalently the orbital period). In order to com- 
bine short coherent segments over increasingly large separa- 
tions in time, one must have an increasingly accurate determi- 
nation of the orbital phase. In short, the longer the observation 
span the more accurately one can determine the orbital period 
and hence more templates must be used to cover the corre- 
sponding search dimension O. 



C. The coherent and semi-coherent metric for long (AT » P) 
coherent segments 

If the uncertainties in each of the parameter space di- 
mensions is sufficiently small or the available computational 
power is large enough then one may be able to work in the 
regime Q.AT » 1 . In this case one would be able to construct 
a search for which it is possible to achieve a coherent obser- 
vation length greater than the length of the orbital period. The 
metric components can then be approximated by taking the 
limit in the regime Q.AT » 1 (given in Appendix |C|l. Note that 
in this regime the physical coordinates 6 - {v, a,f asc ,Q, k, 77} 
(as used for the semi-coherent metric) are a sensible choice 
upon which to base our metric. The metric on a coherent data 
segment can then be well approximated by 



2AT 2 



12v 2 
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(30) 
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where non-zero off-diagonal terms have been approximated as 
equal to zero. We justify this in the same way as done for our 



similar approximation to the semi-coherent metric in Eq. 27 



Just as in the semi-coherent case, the off-diagonal terms con- 
stitute negligible correlations and their inclusion leaves the 
metric determinant essentially unchanged. 



In Fig. [4] we see the degree of agreement between mis- 
matches predicted by the approximate metric and those gen- 
erated by signal injections (using Eq. [2] as the signal model). 
Again we see that our diagonal metric approximation is val- 
idated by the clear lack of correlation between parameters 
and that the discrepancy between approximated and simulated 
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FIG. 2: The results of a simulation in which the mismatch /j(0, A.0) has been computed over a region in the 6-dimensional parameter space 
around the location of a simulated signal. The simulation modeled 500 individual coherently analyzed segments each of 200 seconds length 
spread over a period of t = 10 s seconds. The signal parameters were v = 200, a = 1, Q = 10~ 4 rads s -1 , e = 0.02 and uj = 1 (f asc chosen 
randomly from the range [-P/2, P/2]). Each panel shows the mismatch as a function of each pair of parameters where the mismatch has been 
maximized over all other parameters. The solid black and grey contours indicate the measured mismatch ranging from \i = 0.1 (black) in steps 
of 0.1. The dashed red ellipses are those calculated using the approximate metric given in Eq.|27| 



mismatches in this case is <0.1%. Note that unlike the metric 
on the u parameters this metric is independent of the epoch 
of each segment and hence a template bank generated from 
Eq.|30]can always be identical in equal length segments. This 
type of implementation, as we will see in the following sec- 
tion, can make the combination of results from many segments 
a far easier task than for a search using the u coordinates. In 
addition, it makes the construction of the semi-coherent met- 
ric very simple. Since the semi -coherent metric is equal to the 
"average" metric (averaged over segments) and in this case 
the coherent metric remains unchanged for equal length seg- 



ments, the semi-coherent metric is 



G(0) = lim g(ff) 

AT»P 



and is therefore identical to the coherent metric given in 

We can again compute the scalings associated with the 
number of templates required to cover our 6-dimensional pa- 
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FIG. 3: An orbital parameter sub-space represented as the eccentric- 
ity versus the product of the signal frequency and projected semi- 
major axis. The dashed lines correspond to the semi-coherent met- 
ric's requirement for templates in the eccentricity parameters k and 
j]. Each line represents a different coherent observation length mea- 
sured in units of SIAT. Above each dashed curve a semi-coherent 
analysis would require templates in the eccentricity parameters for 
a mismatch /i = 0.1. Correspondingly, as the coherent integration 
time increases (relative to the orbital period) a limit is achieved at the 
boundary of the shaded region. Below this limit one need not include 
eccentricity in the signal model. 



rameter space as done in Eqs. 28 and 24 giving us 
9V2 



N cn = Me 
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(32) 

where we note that as done in the previous template num- 
ber estimates, we assume that each dimension has a thickness 
greater than the width of a single template span in the corre- 
sponding dimension. For the circular orbit case we obtain 
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In both eccentric and circular orbit cases we have assumed a 
sub-optimal hyper-cubic lattice of templates. As stated in the 
previous section, reductions of factors of a;6.8 and ^2.8 can 
be achieved in the numbers of templates for the eccentric and 
circular orbit cases respectively if an A* lattice is used. 

We note that the number of templates in both the eccentric 
and circular orbit cases does not depend on the total obser- 
vation span t. This means that there is no "refinement" in the 
template bank as one uses larger numbers of segments. In con- 
trast to the previously investigated regime AT <K P, if there is 
an offset in the orbital period then the mismatch from one seg- 
ment will be identical to the mismatch from all other segments 
which will consequently be identical to the total mismatch. In 



short, this is because there are no parameters in our model for 
which fractional losses in SNR accumulate beyond a single 
orbital period e.g. if there are offsets in all parameters, the 
fractional SNR loss after one orbit will be identical to the to- 
tal fractional SNR loss after the next. This is in contrast to the 
behavior of the metric and mismatch in the Ar <K P regime 
discussed in Section |Il A| However, for the same volume of 
parameter space, a significantly larger number of templates 
would be required in the coherent stage for segment lengths 
that far exceed the orbital period compared to the case when 
the orbital period far exceeds the segment length. 



III. DISCUSSION 

We have presented a brief review of the phase model for 
continuously emitting sources in binary systems together with 
a review of the parameter space metric defined on the frac- 
tional loss in SNR when offset from the true signal location in 
parameter space. We then proceeded to apply the latter to the 
former to construct the semi-coherent parameter space met- 
ric for application to searches for such systems in noisy data. 
In addition we provide the outline for the practical applica- 
tion of a search strategy employing the semi-coherent metric 
in the opposing limits that the segments that are analyzed co- 
herently span small fractions of the source orbital period and 
have durations far greater than the orbital period. 

We have shown that in the short duration segment regime 
the semi-coherent metric on the physical parameters 6 be- 
comes approximately diagonal and hence the correlations be- 
tween parameters are removed. This we attribute to the aver- 
aging procedure applied to the metric elements and the spe- 
cific behavior of these elements as a function of orbital phase 
corresponding to each segment epoch. Whilst the metric it- 
self is still a function of location within the parameter space 
making lattice based template placement difficult, the lack 
of parameter correlations does somewhat simplify the proce- 
dure. We note that for the short duration coherent segment 
scenario there is only one term in the semi-coherent metric 
that depends on the semi-coherent observation span. This 
is the element that determines the precision with which the 
orbital angular frequency must be sampled (and conversely 
the precision with which it can be determined in the limit 
of low signal-to-noise ratio). For all other parameters the 
equivalent precisions depend only upon the coherent obser- 
vation time, the signal frequency, orbital angular frequency 
and semi-major axis. This result is identical in nature to that 
obtained in [19| where it was found that only one parame- 
ter in their isolated pulsar model, the pulsar spin derivative, 
was responsible for refinement of the template bank. The ex- 
planation is similar for both cases, both the spin derivative 
and orbital period, in each search respectively, determine how 
results from individual segments are combined whereas the 
other parameters do not. It is only the parameter Q that is ca- 
pable of causing an ever increasing offset in the orbital phase 
(just as the spin derivative in |[T9l is capable of causing a sim- 
ilarly increasing offset in signal phase). As such, coherent 
templates within increasingly separated segments would ex- 



12 




FIG. 4: The results of a simulation in which the mismatch fi(A) has been computed over a region in the 6-dimensional parameter around the 
location of a simulated signal. The simulation modeled a single coherently analyzed segment of 10 6 seconds length containing a signal with 
parameters v = 200 Hz, a = 1 sea, Q = 10 -4 rads sec -1 , sea, e = 0.005 and a> = 1 rad (r asc chosen randomly from the range [-P/2, P/2]). Each 
panel shows the mismatch as a function of each pair of parameters where the mismatch has been maximized over all other parameters (these 
are not slices). The solid black and grey contours indicate the measured mismatch ranging from /j = 0. 1 (black) in steps of 0.1. The dashed 
red ellipses are those calculated using the approximate metric given in Eq.|30| 



perience increasing mismatch for a constant offset in Q. In 
contrast, an offset in the time of orbital ascension (assuming 
an exactly known orbital period) would cause the orbital phase 
to be offset by a constant amount in all segments. Conse- 
quently the mismatch would also be constant and hence the 
average mismatch, equivalent to the total semi-coherent mis- 
match, would also be constant with respect to the observation 
span t. The consequence of this refinement on the orbital an- 
gular frequency parameter is that the cost of the semi-coherent 
stage of the analysis can dominate the cost of the coherent 
stage. However, one is always able to tune the relative costs 
by choosing a shorter or longer segment length. We leave 
the analysis of how to optimally make such choices for future 



work. 

When the duration of the coherent segments are much 
greater than the orbital period there is no dependence upon the 
total observation span and the semi-coherent metric is simply 
equal to the coherent metric since it is constant for each seg- 
ment. In addition, just as for the short segment regime we find 
that the metric is diagonal and that there are therefore no pa- 
rameter correlations. The independence of the total observa- 
tion span on the metric can be explained by the fact that there 
is no dependence upon any of the parameters in the process of 
combining the results from each segment. In contrast to the 
short segment regime, one need not know the orbital angular 
frequency to map a set of physical parameters in one segment 
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to the corresponding (identical) set in another segment. This 
does not mean that parameter estimation is not enhanced with 
the addition of more segments over a longer observation span. 
In this case it is the semi-coherent accumulation of SNR that 
improves parameter estimation whilst the template density re- 
mains constant. 

One interesting point to make given this result is the com- 
parison between the different template scalings with the ob- 
servation parameters for the short and long coherent dura- 
tion regimes. For a fixed set of boundaries in the full 6- 
dimensional parameter space there appears to exist an obser- 
vation span at which it becomes favorable (in terms of using 
fewer semi -coherent templates) to expand the coherent obser- 
vation time to be consistent with the long coherent regime. 
This is because the template scaling with AT is steeper for the 
short segment compared to the long segment scenario. The 
point at which this occurs is dependent upon many interde- 
pendent factors including the volume of the parameter space, 
the relative uncertainties in each dimension and how many 
orders are included in the Taylor expansion in the u coordi- 
nates. This statement is also under the assumption that the 
majority of the computational cost goes into computing the 
semi-coherent detection statistic and that the cost of computa- 
tion of the individual segment statistics is negligible and that 
the computational cost of a long duration template is equal to 
that of a short duration template . It also neglects the issue 
that for a short coherent duration some dimensions are likely 
to not require templates. This occurs when the template span 
is greater than the parameter space range in that dimension. 
Note for example the metric elements corresponding to the 
eccentricity parameters k and rj for the short and long segment 
regimes. The template spacing in a given dimension is propor- 
tional to the inverse of the square-root of the corresponding 
metric element. The ratio between the short and long segment 
k template spacing is V3/(OAr s ) where AT S here is the short 
segment duration. Dependent upon the uncertainty in k, this 
may mean that templates are simply not required in the k di- 
mension for the short segment but are required in the long seg- 
ment regime. If we look at all the contrasting template sizes in 
the short and long duration coherent segments regimes we see 
that in all dimensions except the spin frequency v and the or- 
bital angular frequency Q., a template span is greater by a fac- 
tor oc 1 /(Q.AT) » 1 . The same factor for the spin frequency 
metric components is simply equal to the ratio of the short and 
long observation lengths Ar s /Ari and for templates in the or- 
bital angular frequency dimension the same ratio is given by 
Qt(AT s / AT\. This final result gives the false impression that 
by using short coherent segments one is able to better deter- 
mine the orbital angular frequency of a source than by using 
long segments. We stress here that template placement is re- 
lated to, but not equivalent to, parameter estimation. Template 



placement using the metric defined on fractional SNR (oc the 
log-likelihood) loss determines how finely a parameter space 
must be sampled in order to avoid missing a potential signal. 
Parameter estimation makes use of the likelihood itself and 
it's accuracy improves with SNR. 

We mention at this point that the reader may question our 
choice of phase model parameterization. It is true that pa- 
rameterizations like x = acos(f2f asc ),y = asin(Qf asc ) would 
potentially simplify the metric calculations and that a rescal- 
ing of the frequency parameter could do the same. However, 
since a general "fix-all" re-parameterization eludes the author, 
it was decided that the physical parameters were justifiable 
choices. The fact that the metric elements are functions of the 
parameters themselves is problematic from an implementa- 
tion sense since this will potentially make template placement 
more complicated. Conversely, the parameter choice leads 
to a clearer understanding of how the semi-coherent metric 
scales with the unambiguous physical parameters of the sys- 
tem. 

In conclusion, in deriving the semi-coherent metric for 
known binary systems we have provided a necessary tool for 
the continued sensitivity improvement of searches for contin- 
uously emitting sources of EM and GW emission. With fo- 
cus on GW searches the advanced detector [34] era is fast 
approaching and LMXBs are an ever more realistic detec- 
tion target. The optimality of search sensitivity at fixed com- 
putational cost has been approximately achieved by the Ein- 
stein@Home GW searches [35 1 for unknown isolated neutron 
stars (NSs). The achievement has arisen through the appli- 
cation of semi-coherent search strategies described in |fT9ll 
and are similar to those outlined in this work. It is our hope 
that search strategies based on the semi-coherent methods de- 
scribed here will be employed in searches for the prime GW 
LMXB targets such as Sco X-l and Cyg X-2. 
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Appendix A: Re-parameterized phase coordinates 



In Section II A there is reference to a re-parameterization of 
the phase a ppro ximation given by Eq. |4]into the u coordinates 
used in Eq. 20 The « (m) coordinates for the ra'th segment are 
mapped onto the physical parameters 6 via 
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~(m) _ *0 , („) , 
°0 - ^ + V ( f + f mid - f rf) 



-ov (| (3 + cos[2(f - f asc + f^)Q]) - sin[(f - f asc + )Q] - ^ sin[2(f - f asc + , (Al) 

M ( j m) = y + fl vQ(cos[(f - f asc + f^ d )0] + *cos[2(f - f asc + f^ d )H] + 17 sin[2(f - f asc + ^ d )Q]), (A2) 

4 m) = avQ 2 (27? cos[2(f - f asc + f^ d )Q] - sin[(f - f asc + f^ d )Q] - 2k sin[2(f - f asc + ^ d )0]) , (A3) 

w< m) = -avn 3 (cos[(f - f asc + t^JQ] + Ak cos[2(f - f asc + f^ d )Q] + 477 sin[2(f - f asc + , (A4) 

u ( f = -avQ 4 (81/ cos[2(f - f asc + f« )£2] - sin[(f - f asc + f^ d )Q] - 8* sin[2(f - f asc + )£2]) . (A5) 



Note that we give expressions only up to 4 th order in time. If {v, a, f asc , Q., k, 77} in the limit Ar » P and Q.AT » 1 are 

higher orders are required then it is likely that the coherent 

observation time spans a large fraction (or more) of a single 

orbital period. In this regime using the Taylor-expanded phase 

may not be a sensible choice and it may be more practical to 

use the fully coherent metric on the physical parameters. 

Appendix B: Coherent metric elements for long (Ar » P) 
segment lengths 



To leading order in k and rj (where appropriate) the elements 
of the coherent metric on the physical parameters space = 
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Appendix C: Semi-coherent metric elements the set of parameters 6 = {v, a,f asc ,£2, k, tj}. The following 

elements are valid in the limits of e <k 1 and for t » AT and 
Here we give the expressions for the metric elements in the t » P, 
eccentric orbit case for a phase model defined by Eq. [4] on 
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The following metric elements are the Taylor-expanded versions, to leading order in Q.AT and in e, of those given above 
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The metric given in Eq. 27 is a further approximation in which 



the off-diagonal terms have been set to zero. This approxima- 
tion is valid due to the fact that these terms cause only minor 
mismatch ellipse rotations in the physical parameter space. It 



follows that the diagonal elements themselves are a very good 
approximation to the eigenvalues of the metric and the eigen- 
vectors are very well aligned with the physical dimensions of 
the parameter space. 
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